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ABSTRACT 

We present a generalised formalism for treating the porosity-associated reduction 
in continuum opacity that occurs when individual clumps in a stochastic medium 
become optically thick. As in previous work, we concentrate on developing bridging 
laws between the limits of optically thin and thick clumps. We consider geometries 
resulting in either isotropic or anisotropic effective opacity, and, in addition to an 
idealised model in which all clumps have the same local overdensity and scale, we 
also treat an ensemble of clumps with optical depths set by Markovian statistics. 
This formalism is then applied to the specific case of bound-free absorption of X- 
rays in hot star winds, a process not directly affected by clumping in the optically 
thin limit. We find that the Markov model gives surprisingly similar results to those 
found previously for the single clump model, suggesting that porous opacity is not very 
sensitive to details of the assumed clump distribution function. Further, an anisotropic 
effective opacity favours escape of X-rays emitted in the tangential direction (the 
'Venetian blind' effect), resulting in a 'bump' of higher flux close to line centre as 
compared to profiles computed from isotropic porosity models. We demonstrate how 
this characteristic line shape may be used to diagnose the clump geometry, and we 
confirm previous results that for optically thick clumping to significantly infiuence 
X-ray line profiles, very large porosity lengths, defined as the mean free path between 
clumps, are required. Moreover, we present the first X-ray line profiles computed 
directly from line-driven instability simulations using a 3-D patch method, and find 
that porosity effects from such models also are very small. This further supports the 
view that porosity has, at most, a marginal effect on X-ray line diagnostics in O stars, 
and therefore that these diagnostics do indeed provide a good 'clumping insensitive' 
method for deriving O star mass-loss rates. 

Key words: stars: early-type - stars: mass-loss - stars: winds, outflows - radiative 
transfer - line: profiles - X-rays: stars 



1 INTRODUCTION 

Over the past years, it has become clear that in princi- 
ple all standard spectral mass-loss diagnostics of O stars 
are affected by wind clumping, i.e. by the small-scale 
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wind inhomogeneities that should arise naturally from a 
strong, intrinsic instability associated with the radiative 
line -driving of these win ds (the line-driven instability, LDI, 
e.g. lOwocki et al.|[l988l ). If neglected, such wind clumping 
causes standard diagnostics such as Hq, and IR/radio free- 
free emission, which have opacities that depend on the lo- 
cal wind density squared, to overestimate mass-loss rates 
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(for summaries, see iPuls et alJboOSl : iHamann et al-lbOOSl lFI 
The analysis here examines the degree to which X-ray line 
profiles can provide a mass-loss diagnostic that is relatively 
insensitive to clumping. 

In single O stars without strong magnetic fields, X- 
rays are believed to origi nate in embedded wi nd shocks 
associated with the LDI (jPeldmeier et al.l If 9971 ). and the 
broad emission lines revealed by high-res olution X-ray spec- 
troscopy support this basic scenario (|Kahn et all l200fl : 
ICassinelh et al.ll200fl : fCohen et al.ll2006l ). These X-ray lines 
are often observed as blue-shifted and asymmetric, charac- 
teristics stemming fro m attenuation by bound -free absorp- 
tion in the bulk wind (|MacFarlane et al]|f99fl ). As seen by 
an observer, X-ray photons emitted in the receding part of 
the wind travel farther before escape, and thus undergo more 
absorption, than those emitted in the advancing part. 

For optically thin clumps, the amount of bound-free 
absorption is proportional to the local density, and may 
thereby be used to put additional constraints on mass -loss 
rates. Whilst initial analyses (e.g.. [Kramer et al.ll20oj ) re- 
quired very low mass-loss rates to reproduce the observed 
X-ray lines, more recent investigations with a bette r accou nt 
of th e wind opacity and line blends (|Cohen et al.l lioiOal lbl. 
[20l3) show that rates inferred from X-ray lines are consis- 
tent with those derived from other diagnostics, if clumping 
is adequately accounted for in the other diagnostics. How- 
ever, a possible shortcoming of these X-ray analyses is the 
assumption that clumps are optically thin, which if not met 
would lead to an overestimate of the wind opacity, due to 
the principal effect of porosity. 

Wind porosity models aiming to calculate X-ray lin e 
profiles have be e n dev eloped by, e.g., Feldmeier et al.l l|2003h , 
lOskinova erall (|2004) . and lOwocki fc CohenI (|2006l ') (here- 
after OC06). The first two of these studies assumed the 
clumps to be radially oriented, geometrically thin shell frag- 
ments ('pancakes'), leading to a distinct anisotropic form of 
the effective opacity. In OC06, on the other hand, the clumps 
were assumed isotropic to impinging radiation. Whereas ge- 
ometrically thin shell structures are indeed seen in one- 
dimensional (f-D) LDI simulations, first attempts to con- 
struct 2-D LDI models suggest that these shells break up via 
Rayleigh- Taylor or thin shell insta bilities into small clumps 
of sim ilar angular and radial scales (|Dessart fc Owockill2003l . 
[iooH^). But these initial two-dimensional simulations do not 
yet properly treat the lateral radiation transport that might 
couple material, and so the degree of anisotropy of insta- 
bility generated structure in a fully consistent 3-D model is 
still uncertain. 

F rom the diagnostic side, 0C06 (see also 'Co hen et al] 
l2008f ) argued that for porosity to significantly affect the X- 
ray line formation, unrealistically large porosity lengths, de- 
fined as the mean free path between clumps (see Sect. |3}, 
must be invoked. This view is also supported by the above- 
noted recent attempts to derive mass-loss rates directly from 



^ Diagnostic ultra-violet (UV) resonance lines, which have opac- 
ities that depend linearly on density, are directly affected by 
clumping only if individual clumps are optically thick. However, 
recen t results indicate tha t clumps are indeed thick in these 
lines jPrinia fc Maasall2010h . which then can lead to reduced line 
strengths an d underestimates of mass-loss rates if negl ected in 
the analysis jOskinova et al.ll2007l : ISundqvist et al.ll201ll'l . 
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Figure 1. Schematic of photon of direction n impacting a spher- 
ical clump and a radially oriented (r/|r|) shell fragment. The for- 
mer gives an isotropic effective opacity whereas the latter gives 
an anisotropic effective opacity, since the projected surface area 
depends on direction cosine ^ = cos (9 = n ■ r/|r|. 



X-ray diagnostics. On the other hand. lOskinova et al.l ()2006[ l 
have argued that anisotropic clumps enhance porosity ef- 
fects, and lead to more-symmetric line profiles than if as- 
suming isotropic clumps, in general agreement with X-ray 
observations. 

This paper and its sequel (Leutenegger et al. 2011, 
hereafter Paper II) furt her examine thes e issues. Building 
upon previous works bv lOwocki fc CohenI |200I l and OC06, 
Sects. [2] and [3] develop a generalised formalism for synthe- 
sising X-ray lines, including porosity as caused by either 
isotropic (spherical, or randomly oriented) or anisotropic 
(fiattened, radially oriented) clumps. In addition, we gener- 
alise our models to treat an ensemble of clumps of some dis- 
tribution in optical depth, rather than retaining the assump- 
tion that all clumps are locally identical. Sect.|4]then system- 
atically examines synthetic X-ray line profiles and analyses 
porosity effects for isotropic and anisotropic clumps, as well 
as for uniform and exponential clump distributions. We dis- 
cuss how the shape of the clumps affects the line profiles 
in cases where porosity is important, and how this may be 
used to put empirical constraints on the wind's clump geom- 
etry (leaving detailed confrontation with observed spectra 
to Paper II). Sect. [5] presents first X-ray line profiles cal- 
culated directly from LDI simulations, using the 3- D patch 
method first developed bv lDessart fc Owockil ()2002h . Sect. [6] 
gives a physical interpretation of the analytic porosity mod- 
els presented, showing they can be reconciled with a general 
statistical model derived for stochastic transport in a two- 
component Markovian mixture of immiscible fiuids. Finally, 
in Sect. [7] we discuss our results, compare them to other 
studies, and give our conclusions. 



2 OPACITIES IN A CLUMPED HOT STAR 
WIND 

In our phenomenological model, we assume that the opac- 
ities in the bulk wind can be described using a two- 
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Figure 2. Back-lit rendition of randomly generated spherical (upper row) and radially compressed (lower row) clumps in a spherically 
expanding structured wind (columns 2-7), compared to a back-lit smooth wind (column 1) with t* = 1 and the same total mass between 
an onset radius 1.05 -R* and a maximum visualisation radius 5-R*. The terminal porosity length haa/R* increases from left to right, as 
given by the left header, but clump diameters all scale as del = {0-2/y/TT)r. The right header gives the total number of clumps. The white 
circle represents the star, which also radiates with the same surface brightness as the background source. Appendix A gives details on 
how these visualisations were generated. 



component medium consisting of overdense 'clumps' (de- 
noted with cl) and a rarefied 'inter-clump medium'. The 
distribution of X-ray emitters in the shock-heated wind is 
described in Sect. 13.41 We neglect the inter-clump medium's 
contribution to the opacities, an assumption well justified for 
absorption of X-rays in 0-star winds, due to the generally 
low X - ray o ptical depths found for such stars (| Cohen et al.l 
l2010al . l201ll ). The volume filling fraction of the dense gas is 
fv, thus the local mean density is (p) — fvpci, where the 
angle brackets denote spatial averaging. 



2.1 Optically thin clumps 

X-rays emitted in the wind are attenuated by bound-free 
absorption depending linearly on density. The local atomic 
mean volume opacity per unit length is then (x) = '^(p)i 
with mass absorption coefficient k. By requiring that the 
mean density of the clumped wind be equal to the density 
of a corresponding smooth wind model, i.e. that (p) = psm, 
one immediately recognises the well-known result that for 
atomic processes depending linearly on density, the opaci- 
ties are not directly affected by clumping as long as the in- 
dividual clumps remain optically thin. However, also linear- 
density opacities can via a modified wind ionization balance 
be indirectly affected by optically thin clumping, but mod- 
elling such ionisation equilibria is not a focus of the present 
paper. 

As a comparison, for processes depending on the square 
of the density (e.g. He in hot star winds), the opacities are 
always enhanced compared to smooth models, by a factor 
given by the so-called clumping factor /d = (p'^)/{p)^ = 
fy^, where the latter equality holds when the inter-clump 
medium is neglected, whereby one obtain^ (x) oc (p^) — 



^ In this context, we note that the porosity formalism presented 
in Sect. |3] although developed there for the specific case of X- 
ray line attenuation, is applicable also for continuum processes 
depending on (p^ ) , for example thermal free- free emission, simply 
by exchanging the expression for (x). 



2.2 Porosity 

It is important to realise that (x) may be accurately used 
in radiative transfer models only in the limit of optically 
thin clumps, Td <^ 1. If this condition is not satisfied for the 
investigated process, the radiative transfer becomes more 
complex. For continuum processes such as the attenuation 
of X-rays considered here, optically thick clumps lead to a 
local self-shielding of opacity within the clumps, which in 
turn allows for increased escape of radiation through porous 
channels in between the clump^. The essential effect of 
such porosity is that the 'effective' opacity of the medium 
becomes lower than predicted by an optically thin clump 
model. Thus porosity can mimic the symmetrising effects of 
reduced mass loss on the X-ray line profiles. The purpose 
of this paper is to present a formalism for quantifying this 
reduction and to show how clump geometry and distribu- 
tion affect X-ray line profile morphology in the presence of 
porosity. 



3 A POROSITY FORMALISM FOR 
EFFECTIVE OPACITY 

In analogy with the atomic opacity, we may write the effec- 
tive opacity per unit length of a clump ensemble as (e.g., 
iFeldmeier et al.ll2003l ) 

Xeff = UclAclP, (1) 

where n^i is the number density of clumps, ^ci is the pro- 
jected area (the geometric cross-section) of a clump for di- 
rection n, and P is the probability that a photon impacting 
a clump gets absorbed. This probability obviously depends 

^ For line formation in a rapidly accelerating clumped medium, 
optically thick clumps lead to corresponding velocity gaps, 
through which lin e photons may escape without ever i nteracting 
with the material dOwockilbOol : ISundgvist et al.ll2010l) . This is a 
consequence of the Doppler shift, leading to a picture wherein the 
clump length scales can be comparable to (or even larger than) 
the extent of the lines' resonance zones, which is a limit wherein 
the porosity formalism developed here is not directly applicable. 
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on the clump optical depth, P — P{tci), which we here 
characterise by the average over all possible chord lengths I 
through the clump for rays of direction n, 



Ac 



A, 



Ac 



(2) 



where Md is the mass of the clump and the last equality 
assumes that k is constant over the clump. 

Another useful quantity is the local mean free path 
between clumps, a lso known as the porosity length h 
lOwocki et al.ll2004l '). 



h ; 



riciAci 

Using this definition, Tci may be written 
kMci 



Tcl = 



Aci 
whereby 

P{ra 



= n{nc\Mc\)h = K{p)h = (x)^. 



Xcff 



h 



Xoff 

(x) 



Tcl 



(3) 



(4) 



(5) 



For P = 1, Eq. [S] returns the atomic opacity independent 
result Xeff — 1/h, demonstrating how the porosity length h 
can also be interpreted as a photon's mean free path in the 
limit of only optically thick clumps. Further, for a constant 
clump density Eq. [2]yields Tci = Kpdlav, with average chord 
length £av, which then recovers the commonly used form 
h = Isuv/ fv for the porosity length. 



3.1 Isotropy vs. anisotropy 

Note that the porosity length as defined in Eq. |3]is a strictly 
local quantity. And because Td = {x}^^ this means that the 
(an)isotropy of the effective opacity (Eq.[5]) depends only on 
Aci, and is independent of the spatial variation of h associ- 
ated with the global wind expansion. Thus, spherical clumps 
(Fig. [1]), as well as randomly oriented clumps of arbitrary 
shape, will have an isotropic effective opacity. For these cases 
then, assuming spherical symmetry for the global wind ex- 
pansion, one may set h — h(r) for all directions n impacting 
the clump. 

However, now let us consider a specific wind model 
in which the clumps really are randomly distributed, but 
radially oriented, geometrically thin shell fragmen ts ('pan- 
cakes') ijFeldmeier et al.ll2003l : IOskinova et al.ll2004l ). In such 
a model, the projected clump area is Acin- v/\r\ = A^^il/ij, 
where we here identify the projected area for a radially di- 
rected ray with the area for an isotropic clump (Fig. [1} . 

This implies the effective opacity retains its basic form 
(Eq. [5} also for such anisotropic clumps, but the clump op- 
tical depth becomes larger for oblique rays. 



T-A{r,p) = {x)h{r,p) = {x)h{r)/\p\, 



(6) 



where the isotropic case is recovered by setting p = 1, 
and the radial dependencies of the mean opacity and di- 
rectional cosine have been suppressed. This paper considers 
only isotropic and 'pancake' geometries, but the porosity 
formalism outlined above applies generally to any clump ge- 
ometry described by Aci = ./(n), where / is some function. 

To illustrate such isotropic vs. anisotropic absorbing 
media. Fig. [2] compares a random distribution of spherical 



clumps and radially oriented pancakes, as illuminated by 
a uniform background source (see Appendix A). Note that 
for visual clarity, we extend these visualisations only to an 
outer radius 5 i?* . We stress that Fig. [2] is for general il- 
lustration purposes; the uniform background illumination is 
not a distribution of X-ray emitters. Also, we repeat that 
though assumed in the visualisations, spherical clumps are 
actually not a necessity for obtaining isotropic porosity; the 
requirement for this is rather that the clumps be randomly 
oriented, see above. 

3.2 Bridging laws for the effective opacity 

For attenuation with a given local clump optical depth rd = 
To, let us assume the probability of absorption simply takes 
the basic form P = 1 — e~'^° . While formally exact only in 
cases where all chord lengths across the clump are equal, 
this expression is a suitable approximation that yields with 
Eq.[5] 



To 



(7) 



This 'single clump' bridging law now has the correct values 
in the limiting cases; it returns the atomic mean opacity 
when To <gil and is independent of it when ro 3> 1. And as 
discussed in the preceding section, the bridging law equation 
applies for both isotropic and anisotropic effective opacity 
models, however with different expressions for the clump 
optical depth. 

Assuming the effective mean free path scales as Xc« ~ 
(x)^^ + h, an even simpler 'inverse' (or 'Rosseland', see 
OC06) bridging law was invoked in OC06, 



(X> 1 + TO ' 



(8) 



which also has the correct optically thin and thick limits. 
By considering the direction dependent mean free path, this 
inverse bridging law can be realised also for anisotropic mod- 
els. 

In OC06, the practical motivation for invoking Eq. [8] 
was because the optical depth integral for X-ray attenuation 
could then be solved analytically. However, Sect. [S] shows 
it also happens to represent the bridging law that follows 
from assuming the local clump optical depth distribution 
function obeys Markovian statistics. Thus, the single-clump 
and inverse bridging laws (Eq.Qand Eq.|8} differ in that the 
former assumes all clumps have the same local optical depth 
(for a given direction), whereas the latter averages over an 
exponential distribution in rd. 

3.3 Velocity stretch porosity 

OC06 assumed that the porosity length scales with the lo- 
cal radius, but for mass-conserving clumps, such a scaling 
is only appropriate for isotropic expansion. For clumps re- 
leased into a radially expanding stellar wind, the wind ac- 
celeration will 'str etch' the clump spacing in proportion to 
the wind velocity (|Feldmeier et a l. 2003). The analysis here 
assumes this velocity stretch form for both isotropic and 
anisotropic porosity. 

This distinction is most easily seen for the radially frag- 
mented shell model, in which the average radial separation 
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between two shells is Ar — h(r). For shells moving radi- 
ally according to a '/?- velocity law', w{r) = v{r)/voo = 
(1 — Ri,/r)^ , where 7?* is the stellar radius, the separa- 
tion is h(r) = haowir), with the parameter h^o representing 
the asymptotic radial separation as w — >■ 1. For simplic- 
ity, this paper assumes the prototypical value 13 = 1. The 
quantity /ioo/«oo represents the average time between two 
consecutive shell pa ssings at a fixed radial point in the wind 
l|Sundavist et alJlioiQ ) , which may also be in terpreted as the 
inver s e of a 'fragmentation frequency' no (|Oskinova et alj 
l2oo4 booel v For example, hoo = R* gives a fragmentation 
frequency no = Vaa/R* that is equal to the inverse of the 
wind flow time. 



3.4 X-ray line transfer in porosity models 

To compute X-ray emission line profiles, we solve the stan- 
dard formal integral of radiative transfer, using a customary 
(p, z) coordi nate system and followi ng the basic procedure 
described in lOwocki fc CohenI (|200ll ) for the distribution of 
X-ray emitters. Since our primary interest here is the shapes 
of the lines, all resulting fiux profiles have been normalised 
to a unit maximum, = J^x/Max(J-'x). For simplicity 

we assume that the X-ray emission begins at a certain onset 
radius 7?o = 1.5 7?*, and is constant beyond it. This onset 
radius is consistent with that typically pred icted by con- 
serva tive, self-excited LDI simulations (Runacres fc Owockil 
120021 '). and both 7?o ~ 1.57?* and a constant X-ray filling fac- 
tor are supported by obs ervations (jLeutenegger et al ] |2006l : 
ICohen et atlbOOd I2OI1I ). Moreover, both LDI simulations 
and the observed X-ray luminosities indicate that only a 
very small mass fraction, less than 1 %, of the stellar wind is 
shock-heated to X-ray emitting temperatures at any given 
time. Simulations and observed lack of X-ray variability fur- 
ther indicate that there are numerous sites of X-ray emis- 
sion distributed throughout the wind, justifying our assump- 
tion of a smoothly distributed X-ray emitting plasma above 
the onset radius. Note though, that 7?o does not necessar- 
ily equal the clump onset radius 7Zci, which observations 
typically indicate is located much cl oser to the photosphere 
(Puls et al.ll2006l : ICohen et al.ll201ll , see also Fig. [2] for a vi- 
sualisation). But in the velocity stretch porosity models, we 
have verified that the exact value of this 7?ci is not important 
for the resulting X-ray line profiles. 

The absorption of X-rays emitted at position Zc along 
a ray with impact parameter p is given by the optical depth 
integral 



(p, Zc) = / Xo«{z,p)dz 



(9) 



where the effective opacity accounts for any porosity. The 
opacity in a smooth or optically thin clump model, due 
purely to the atomic mean opacity (x), is proportional to 
the mass-loss rate M of the star and here characterised by 
a fiducial optical depth r* = Mk,/ (47r7?*i;oo), with wind ter- 
minal speed Voa- To evaluate Eq. [5]for our porosity models, 
the only additional input parameter required is the poros- 
ity length h (see Eqs. I4I8|) . This holds for isotropic as well 



as anisotropic effective opacity models, and for the single- 
clump as well as the inverse bridging laws0 



4 X-RAY LINE PROFILES FROM ANALYTIC 
POROSITY MODELS 

Fig. displays synthetic X-ray line profiles calculated using 
the four possible combinations of isotropic vs. anisotropic 
effective opacity and single-clump vs. inverse bridging laws. 
The figure clearly shows that profiles calculated using the 
two different bridging laws are very similar (see also Fig. [5l 
as well as Fig. 1 in OC06), despite representing two very 
different clump optical depth distributions (Sect. [5^ . This 
indicates that the effects of porosity on X-ray line profiles 
are not very sensitive to the specific local distribution in t^i- 
We discuss this important result further in Sect. [S] 

The second key feature of Figs. [3] and [5] is the promi- 
nent 'bump' visible close to line centre in profiles calculated 
with anisotropic effective opacity. Conceptually, we may un- 
derstand this a s a 'venetian blind' effect (Fig. U see also 
IFeldmeier et all 2003); since the fragmented shells are radi- 
ally oriented, the blinds are closed for radial photons, but 
open up for more tangential ones. This leads to increased 
escape for photons emitted close to line centre, since the 
line emission wavelength scales with direction cosine ^ as 
(e.g.. lOwocki fc CohenllioOll l. 

Another way to look at this effect is to consider the 
optical depth integral for anisotropic effective opacity in the 
Tci 3> 1 limit. 



r(p, 2) 



XcH{z,p)dz : 



h{r. 



■dz, 



(10) 



which shows that, since dr = fj,dz, the optical depth in this 
limit is set simply by counting up the number of porosity 
lengths. In the plane-parallel limit of radially oriented, geo- 
metrically thin but optically thick fragments, all tangential 
= 0) photons would escape. However, due to sphericity ef- 
fects (i.e., that jj, increases as the photon propagates through 
the wind) also photons emitted initially in the tangential 
direction will suffer some absorption (see Fig. |3|. Thus the 
end result is not complete transmission, but a characteristic 
bump stemming from the reduced integrated optical depth 
for photons emitted around 2: ~ 0. For isotropic porosity, on 
the other hand, no fi factor enters in Eq. 1101 and therefore 
no bump appears in these profiles. 

This quite distinct and systematic difference in the 
shape around line centre between models with isotropic and 
anisotropic effective opacity is a key result of the present 
analysis. Indeed, one can us this difference to set empirical 
constraints on the clump geometry by confronting synthetic 
X-ray spectra with observed ones, as will be done in Paper II. 

Generally, Fig. |3] confirms earlier results by OC06 that 
in order to achieve a significant effect on the profiles, rather 
large porosity lengths, hoo > R*, are required. However, 



^ Source codes to all X-ray line porosity models presented 
in this section are publicly available at (package windprof) 
'http: / /heasarc.gsfc. nasa.gov/docs/xanadu/xspec/ models/wind 
prof.html, where models for the broadband absorption of X-rays 
(package WINDTABS, ILeutenegeer et al.ll201Cll) also can be found. 
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Figure 3. Synthetic X-ray line profiles for the 'single clump', Xoft/(x) = (1 ~ ^■ii'i the 'inverse', Xcff/(x) = 1/(1 + effective 

opacity bridging laws and for different porosity length parameters hco/R*, as labelled. All models assume an onset radius for the X-ray 
emission Rq = 1.5iJ*. The abscissae display the dimensionless wavelength x = (A/Aq — l)c/i;oo, with Aq the line-centre wavelength, and 
the ordinates display the normalised flux. Black, red, blue, and green dashed lines have an optical depth parameter Tt = 0.01, 1, 5, 10, 
respectively; for a non-colour separation, an increased r* means a more blue-shifted peak flux. Note that we have set hao = (the 
uppermost panel) to be equivalent to assuming only optically thin clumps. 



Fig. [S] reveals that for a value of = 2.5 , representa- 
tive o f the prototypical O supergiant Pup (| Cohen et al.l 
l2010al ). the anisotropic porosity model displays significantly 
higher flux around line centre than the other profiles, also 
for hao = R*- Further comparisons show that, for the pa- 
rameters used in Fig. O the isotropic porosity model is well 
matched by a corresponding optically thin clump model with 
r* reduced by ~ 20 percent (a detailed quantification of this 
trade-off will be provided in Paper II). In contrast, such a 
simple optical depth reduction does not reproduce the line- 
centre region in anisotropic porosity models, i.e. there is 
no simple trade-off, or degeneracy, between mass-loss rate 
and anisotropic porosity. Thus, whereas it will be difficult 
to distinguish between optically thin clumps and moderate 
isotropic porosity with somewhat higher r^'s, careful line- 
fitting to observations should be able to identify, or refute, 
an anisotropic porosity. 



in a large portion of the wind, the profile shapes do indeed 
become quite independent of atomic opacity. In the bottom 
panel of Fig. |3l profiles computed using an anisotropic ef- 
fective opacity are actually very similar to the flat-topped 
proflles stemming from computations without any X-ray ab- 
sorption, whereas proflles computed using an isotropic ef- 
fective opacity retain a certain degree of asymmetry. This 
result differs from that found in OC06, wherein very near 
symmetry was achieved for isotropic proflles with very large 
porosity lengths, and comes about due to the h (x v scaling 
adopted here, which implies signiflcantly shorter porosity 
lengths in the lower wind regions than the h ot r used in 
OC06. However, for more moderate values of /loo (top three 
panels in Fig. |3}, all proflles are opacity dependent, as ex- 
pected since in the formation of these lines, porosity is a 
secondary effect. 



Finally, for very large porosity lengths, giving Tci ^ 1 
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Figure 4. Illustration of the 'Venetian blind' effect seen in poros- 
ity models using an anisotropic effective opacity. The dashed ar- 
rowed lines represent two different p-rays and the observer is as- 
sumed to be located at Zrx,. 



Optically thin clumping 
Anisotropic, inverse 
Isotropic, inverse 
Isotropic, single clump Z 




Figure 5. Line profiles for hao = R* and r* = 2.5, using different 
effective opacity laws, as labelled. 



CAK 

LDI, 1 degree patch 



— LDI, 3 degree patch 




Figure 6. Line profiles for t* = 2.5, calculated from a smooth 
CAK model and structured LDI models with patch sizes I and 3 
degrees (see text), as labelled. 



5 X-RAY LINE PROFILES FROM LDI 
SIMULATIONS 

The above simplified analytic models demonstrate that the 
effect of optically thick clumps on wind continuum absorp- 
tion can be well characterised in terms of the wind porosity 
length, parametrised here by the asymptotic value hoo and 
assuming a radial variation set by the velocity-stretch form, 
h{r) = ha^{l- Ri,/r) (Sect. lO)) . 

But inspection of LDI simulations suggests a substan- 
tially steeper radial variation for the separation between 
instability-generated wind clumps. Even in 1-D models in 
which the separation can become of the order of a stellar 
radius in the outer wind (implying hoo ~ R*), the initial 
clump structure formed near the onset radius r ~ 1.5-R* 
tends to have a much smaller separation scale, of the order 
of the Sobolev length Lsob ~ {vth/v)R* ~ O.Oli?*, where 
vth is the ion thermal velocity. The sharp increase in separa- 
tion from this onset comes not just from velocity stretching 
from the overall wind acceleration, but also from collisional 
merging of clumps with substantial radial velocity disper- 
sion. For a given asymptotic porosity length /loo, LDI mod- 
els thus tend to have smaller inner-wind porosity lengths 
than assumed in the simple velocity stretch scaling. Since 
it is in this inner wind region that clumps can become op- 
tically thick for X-rays, this suggests that LDI models will 
show even weaker porosity effects than implied by the an- 
alytic profiles shown in Figs. |3]and[5l To demonstrate this 
explicitly, we now present some first sample calculations of 
X-ray profiles computed from 1-D LDI simulations that are 
phased randomly among 3-D patches of a parametrised an- 
gular size (as first developed in .Dessart fc Owocki 200^). 
The details of this patch geometry as implemented in ou r 
radiative transfer code are given in ISundavist (|201ll ). 
Here we adapt this code to synthesise X-ray line profiles by 
making the following assumptions: The X-ray emission is 
assumed to have a fixed spatial form independent of wind 
structure, scaled in proportion to the density squared of a 
smooth CAK wind and with an onset radius Ro = 1.5/?*. 
The bound-free absorption is then calculat ed directly from 
the st ructured LDI simulati on presented inlSundavist et al.l 
(|201lh (computed following .Feldmeier et al.l 19971" 

Fig. |6]plots X-ray line profiles for the same wind optical 
depth used for analytic models in Fig. (5] namely r« = 2.5. 
The curves compare a smooth CAK model to structured 
LDI models with patch sizes of 1 and 3 degrees. The overall 
shapes agree well with corresponding non-porous analytic 
models, except for small differences due to the fact that the 
CAK velocity law does not exactly follow the phenomeno- 
logical /3 = 1 law. But the key point is that the LDI profiles 
are very similar to the CAK profile, implying little porosity 
effect; however, we note that for a 3-degree patch size, there 
is a small, but noticeable, 'bump' around line centre, pre- 
sumably associated with the anisotropic (pancake) nature of 
the clumped structure in such models. 

These results confirm that porosity is as marginal 
in LDI simulations as in corresponding analytic stretch- 
porosity models with comparable asymptotic clump sep- 



arations (i.e., hg 



Ri, ). Mor eover, in 2-D LDI models 



IIDessart fc Owocki|[200^ . l2005al ). clumps can also be bro- 
ken up by shearing and associated effects, leading generally 
to more, smaller, and less optically thick clumps, charac- 
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terised by even smaller porosity lengths. Future work will 
implement our radiative transfer tools also in such genuinely 
multi-D instability models. 



6 A PHYSICAL INTERPRETATION OF THE 
EFFECTIVE OPACITY BRIDGING LAWS 

As already noted, a basic difference between the two bridg- 
ing laws adopted in this paper is that one (Eq.[71 the 'single 
clump' law) assumes a locally constant clump optical depth 
for all clumps, whereas the other (Eq. [8l the 'inverse' law) 
represents a certain distribution in Td. Before discussing 
porosity effects on the synthetic X-ray line profiles in Sect.[7l 
this section examines the nature of this distribution. 



6.1 Distribution laws for 

To investigate how a distribution of clump optical depths 
affects effective opacity scalings, let us assume that the ratio 
of the effective to the mean opacity scales as the ratio of an 
effective clump optical thickness to its mean. 



Toff 

{r) 



-dr 



(11) 



where we for clarity have dropped the indices on the clump 
optical depths. Eq. [11] introduces /(r), the normalised 
distribution function of clumps, and t^s, the distribution 
weighted mean of the clump optical depth. Selecting a 
weighing function (1 — e~'^)/r ensures that the single clump 
bridging law (Eq.O is recovered from Eq. Illl when the dis- 
tribution function is a Dirac delta function, 5(t — tq). 

Let us now choose a specific distribution function of the 
exponential form 



fir) 



TO 



(12) 



which has a mean value tq . Using Eq. [12] in Eq. [TT] yields 
directly the inverse bridging law Eq. [S] but with tq now only 
representing the mean clump optical depth, rather than a 
unique one as in the exponential bridging law Eq. [7] Again, 
despite the large difference between this clump distribution 
and the one assuming a constant r, the two bridging laws 
give similar results, as demonstrated in Sect. [4] 



6.2 Connecting Markovian statistics to 
exponentially distributed clumps 

We now show that the bridging law resulting from this ex- 
ponential distribution turns out to be a special case of a 
general scaling derived for a stochastic mixture of two fluids 
that follow Mark ovian statistics (e.g., Lcvcrmorc ct al. 1988; 
IPomranineil 199lh . The Markov assumption is that the future 
state of the system only depends on its present state, and 
not on its history. Along any given ray through the medium, 
if the fluid is of component at location s, the probability 
of it being of component 1 at s -|- ds is Po,ids, where Po,i 
is independent of how far back along the ray the last tran- 
sition (from fluid 1 to 0) occurred. Under this assumption, 
the length scales traveled within the fluid components are 
random variables described by Poisson distributions, with 



Po,i identifled as the inverse of £o, the mean distance a pho- 
ton travels along the ray in fluid before finding itself in 
component 1. A similar definition applies for £i. 

For su ch a two-component Mar kov model with opacities 
Xo and yi . iLevermore et all (|l988h derived for the effective 
opacity in the pure absorption case 



Xcff 



ix) + XoXi4 
1 + (poXi +PiXo)4 



(13) 



where pi = ^i/(4 + ti), (x) = PoXo + XiPi is Ihe mean 
opacity, and £c = ^o^i/(^o -|- ^i) is the correlation length 
(|Pomraninglll99ll ). Identifying the clumps in our model with 
component 1, and assuming the inter-clump medium to be 
void (xo = 0), we find for this 'clump-(-void Markov model' 



Xcff 
(X) 



1 +PoXi^ 



l + {xWo l + {x}h' 



(14) 



Here po represents the probability that a photon is in the 
void medium, while £o is the distance the photon travels 
in the void before encountering a clump; the product poio 
thus represents the photon mean free path in the case of 
optically thick clumps, which is also the porosity length h, as 
given by the final equality in Eq. 1141 Comparison with Eq. [8] 
then shows that the effective opacity bridging law for an 
exponential clump optical depth distribution is equivalent 
to that for this statistical clump-|-void Markov model. 

Indeed, recalling that the Markov transport is defined 
along a given ray, we may make the same identification for 
the anisotropic porosity model, with the photon mean free 
path along the ray then being scaled by 1/|m|. 



6.3 Exponentially truncated power-law 
distributions 

While there are not many observational constraints on 
the distri bution of clumps in a hot s tar wind (see how- 
ever Lcpi ne fc Moffat! Il999l : iDessart fc O wocki 2005bll , the 
above identification with the Markov model at least places 
our porosity models on a robust and well-known statisti- 
cal ground. The inverse bridging law should therefore be an 
appropriate standard choice for porosity applications such 
as the X-ray line formation considered here, but perhaps 
also for, e.g., porosity moderated continuum driven wind 
models of stars formally excee ding th e Eddington lumi- 
nosity, as investigated by Owoc ki et all l|2004l ). Indeed, al- 
though not explicitly studied in that paper, we note that 
the Markov model represents a special case of the exponen- 
ti ally truncated power- law distribution of clumps considered 
in lOwocki et al.l l|2004l ). namely the one with power index 
Op = 2. Thus, reasonable extensions of the two canonical dis- 
tributions studied here could readi ly be done by u sing s ome 
other power index variant given in lOwocki et al] (|2004l '). 



7 DISCUSSION AND CONCLUSIONS 
7.1 Isotropy or anisotropy 

Let us next compare our analysis to that bv lOskinova et all 
(|2006l) . These authors also pointed out the differences be- 
tween isotropic and anisotropic effective opacity and car- 
ried out a comparison, however only for the specific case 
of — 10 and a fragmentation frequency (see Sect. I3.3|l 
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no = 1.4 X 10"^ s~^. Tak ing the paramet ers for the O su- 
pergiant Pup adopted in lOskinova et ahl . this corresponds 
to hoo = Voo/no = 4.24J?*, which in turn roughly corre- 
sponds to the bottom panel in our Fig. [3] Indeed, the pro- 
fil es displayed in t hat panel agree well with those in Fig. 16 
of lOskinova et aU : both figures illustrate that for such very 
large porosity lengths, profiles computed using anisotropic 
effective opacity are nearly symmetric (Sect. [4]). 

This comparison suggests an overall good agreement 
among the results found by the different groups. But as 
shown in Sect. U anisotropic porosity line profiles, with 
their characteristic 'bump' at line centre, are qualitatively 
different than isotropic porosity or optically thin clump- 
ing profi les. Thus the good s tatistical fits presented for 
Pup by ICohen et al.l l|2010al 'l. without invoking porosity, 
seem some what contradic t ory t o the good visual fits pre- 
sented by Oskinova et al.l (|2006l ). using models with mod- 
erate anisotropic porosity hoo ~ R*- Paper II will further 
examine and quantify these differences between anisotropic 
porosity on the one hand, and isotropic porosity or optically 
thin clumping on the other. 



7.2 Is porosity important for X-ray line mass- loss 
diagnostics? 

The LDI simulations presented in Sect. [5] indicate small 
porosity lengths and negligible porosity effects on X-ray line 
profiles. Such small porosity lengths also have some indirect 
empiric al support. N a mely, t he mass-loss rate derived for 
Pup bv lCohen et al.l (|2010al ). without invoking porosity, is 
only margi nally lower than the upper limit mass-loss rate 
derived bv IPuIs et al.l (|2006| ). by assuming an unclumped 
outermost radio emitting wind, while allowing for clumping 
in the intermediate and lower wind. Because of the trade-off 
between porosity and mass-loss rate then, if porosity lengths 
large enough to significantly affect the X-ray line profiles 
were to be adopted, hoo > R*, the inferred X-ray mass- loss 
rate would be higher than this upper limit. That is, such 
multiwavelength considerations indicate that a significant 
porosity effect on X-ray based mass-loss rates is incompati- 
ble with diagnostic results from other wavebands. 

Overall, we thus conclude that porosity effects on X- 
ray line profiles are likely to be, at most, a marginal effect 
in typical O stars. This is supported also by the low opti- 
cal depths found for C Pup as well as for the even denser 
wind of HD93129A (jCohen et al.ll2010al . [2OI1I ). Since most 
O stars will have characteristic r«'s significantly lower than 
these, porosity effects should be negligible. The upshot is 
that X-ray line analysis may indeed provide the best avail- 
able 'clumping insensitive' diagnostic of O star mass-loss 
rates. 
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APPENDIX A: GENERATION METHOD FOR 
CLUMPED MEDIUM ILLUSTRATIONS 

The clumped medium illustrations in Fig. [5] were generated 
by following the radial expansion of mass-conserving clumps. 
As with the analogous illustration in Fig. 3 of OC06, we 
assume the clump scale I increases in proportion to the local 
radius, 



l{r) = h 



R* 



(Al) 



However, instead of the OC06 assumption of a purely 
isotropic ('Hubble-law') velocity expansion « ~ r, we now 
use a standard P — 1 wind velocity law. For clumps of 



projected area Ad = / 



and local volume density 



rici oc l/{vr^), this gives the associated radial variation of 
the porosity the desired 'velocity-stretch' form. 



h{r) 



1 v(v) 
— -r- = ftoo-^ = hoo{l - R,/r) . 



(A2) 



For specified clump parameters h and hoo, the clump num- 
ber density is thus given by 



Uciir) 



voo Ri 



hAci hooti r'^v{r) 



(A3) 



Note that, unlike the OC06 isotropic expansion model, the 
clump volume filling factor in this velocity-stretch scaling 
is not constant, but varies spatially as /v oc Udl^ oc r/v(r), 
which actually is quite consistent with derived observational 
constraints (e.g.. |Puls et al.ll2006l ). 

The cumulative number of clumps up to a radius r 
above the clump onset radius R^ is 



N{r) = 



47r / Ucir' dr' 



hooH 



Rc 



R, 



In 



R, 



(A4) 



clumps is A'^tot ~ Af(_Rniax). Since Fig. [2] uses fixed param- 
eters L/R, = 0.1, Rci/R, = 1.05, and R^^^/R, = 5, this 
number scales with I//I00, and ranges from A'tot = 1308 
for the largest porosity length hoo/R* = 8 in the right- 
most column, to Ntot ~ 41 881 for the least porous case 
hoo/R* ~ 0.25 in column 2. 

A random set of A'tot clumps with the required statisti- 
cal distribution in radius can now be drawn simply by gen- 
erating a set of pseudo-random numbers TZi over the range 
[0, 1], and inverting the normalised cumulative distribution 
function F{r) = N(r)/Ntot = Ti-i to find the radius r, 



r(7^i) ^R*{1 + ProductLog[exp(Cci + Cmax7^i)]) 

with the constants given by 

Col = Roi/R, - 1 + \n{Roi/R* - 1) 

Cmax = (-Rmax ^ Rcl) / R* 4^ lu ^- 



Roi - R, 



(A5) 

(A6) 
(A7) 



Likewise, we use the assumed statistical spherical symmetry 
and generate the clump angle coordinates in azimuth (j> and 
colatitudinal cosine p through additional pseudo-random 
numbers, 



(A8) 
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In this system it is most convenient to assume the clumps 
are viewed from above the /j. — 1 pole, with thus /i > 
(p < 0) representing clumps in the foreground (background) 
hemisphere. 

With the random set of clump positions in hand, the 
clump sizes are set by Eq. (|Aip . with all cases in Fig. [2] 
using h/R, = 0.1. To give the associated projected area 
— Ac\ — 7vd^i/4, the clump diameters are set to d^i — 
2//0F. 

For the spherical clumps in the upper row of Fig. (2] the 
transparency of each individual projected clump disk area 
is set by exp(— Tcj), where the surface-averaged clump opti- 
cal depth is Tci = Tt,hcx,R*/r^ ■ For the radially compressed 
clumps in the lower row, the associated clump optical depths 
are increased by their projected areas are reduced 

through foreshortening their radial extent by a factor j/ij, 
while keeping their radially perpendicular extent equal to 
the local clump diameter dd- 

Finally, clumps in the back hemisphere (/i < 0) that 
are behind the star (with r -^/l — /i^ < R^,) are simply not 
drawn. This effectively means the clumps directly in front of 
the star appear as if illuminated by a stellar surface bright- 
ness equal to the back illumination of the clumps outside 
the stellar limb. 
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